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ABSTRACT 

We suggest that the irregular structure in Saturn’s B ring arises from the 
formation of shear-free ring-particle assemblies of up to ~ 100 km in radial ex¬ 
tent. The characteristic scale of the irregular structure is set by the competition 
between tidal forces and the yield stress of these assemblies; the required ten¬ 
sile strength of ~ 10® dyn cm“^ is consistent with the sticking forces observed 
in laboratory simulations of frosted ice particles. These assemblies could be the 
nonlinear outcome of a linear instability that occurs in a rotating fluid disk in 
which the shear stress is a decreasing function of the shear. We show that a sim¬ 
ple model of an incompressible, non-Newtonian fluid in shear flow leads to the 
Cahn-Hilliard equation, which is widely used to model the formation of structure 
in binary alloys and other systems. 

Subject headings: planets: rings — celestial mechanics 


1. Introduction 

The Voyager 1 and 2 spacecraft, which flew past Saturn in 1980 and 1981, revolution¬ 
ized our understanding of the Saturn system. One of the remarkable features discovered by 
Voyager was rich radial structure in Saturn’s rings. Although some of this structure, mostly 
in the outer or A ring, is known to arise from density or bending waves generated by discrete 
resonances with the inner satellites, most of the radial structure remains unexplained, par¬ 
ticularly in the main or B ring. Horn & Cuzzi (1996) point out that wavetrains associated 
with known resonances cover less than 1% of the radial extent of the A and B rings. The 
remaining vast majority of the structure is often called “irregular” since the bright and dark 
features show little or no long-range coherence. 

The purpose of this paper is to suggest and examine a novel explanation for irregular 
structure in Saturn’s rings. Ring particles are likely to have (weak) cohesive forces, and 
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therefore can assemble or “freeze” into structures much larger than an individual ring parti¬ 
cle. The size of these structures is limited by tidal forces from Saturn and collisional erosion 
by impacting particles. We suggest that much of the irregular structure may consist of alter¬ 
nating annuli of “solid” and “liquid” ring material, the former consisting of an assembly of 
ring particles frozen into rigid rotation around Saturn, and the latter consisting of individual 
ring particles in differential rotation. Thus we hypothesize that the irregular structure is 
primarily a manifestation of variations in shear rather than surface density. 

Rich spatial structure associated with the co-existence of two phases is a feature of 
many physical processes, including ferromagnetism, spinodal decomposition in alloys, crystal 
growth, chemical reactions, and even traffic flow. We shall exploit some of these analogies 
in §5. One complication in Saturn’s rings is that the physics is non-local, because the tidal 
stresses on a rigidly rotating annulus depend on its total radial extent. We estimate these 
stresses and the corresponding upper limit to the radial extent of rigidly rotating annuli in 
§4.2. 

We begin by reviewing the observations of the irregular structure in §2 and competing 
theoretical explanations in §3. In §4 we review the equations that describe the dynamics of 
the ring fluid. In §5 we describe a simple toy model that illustrates some of the important 
dynamical behavior associated with the irregular structure, and sets the model in the context 
of the theory of phase transitions. 

The suggestion that some areas in planetary rings are locked into solid assemblies is not 
new to this paper. Wisdom & Tremaine (1988) point out that dense rings may have a liquid- 
solid phase transition; they did not find any such transition in their numerical experiments 
on inelastically colliding hard spheres, but this is not surprising since they did not model any 
cohesive forces. Mosqueira & Estrada (2002) suggest that ring-particle shocks or “jams” may 
occur in sectors of converging flow in an eccentric ring; however, our focus is on axisymmetric 
ring-particle structures held together by cohesive forces, while theirs is on non-axisymmetric 
structures held together by self-gravity and ram pressure. 


2. Observations 

Voyager images reveal irregular structure throughout the A and B rings, over a wide 
range of scales ranging from several hundred km to the resolution limit of about 5 km. In 
regions where the ring is sufficiently transparent, the surface-brightness variations track the 
optical-depth variations detected during stellar occultations observed by Voyager (Horn & 
Cuzzi 1996). However, in the main part of the B ring the optical depth is so high that there 
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is little or no signal from the occultation experiments; moreover, the brightness variations 
in this region cannot be explained by optical depth variations, since the solar elevation at 
the time of the encounter was so low (< 8°) that variations in optical depth have negligible 
influence on the surface brightness in reflected light. Thus, other effects such as particle 
albedo or phase function, which may or may not be correlated with optical depth variations, 
must be responsible for most of the irregular structure in the B ring. Comparison of images 
at different phase angles (Cuzzi et al. 1984) suggests that some features are due to variations 
in albedo and others due to variations in phase function. The surface-brightness ratios at 
different phase angles appear to be bimodal, suggesting that there are only two distinct 
phase functions. 

Most of the irregular structure appears to be axisymmetric, at least on scales > 50 km 
and away from major resonances. At smaller scales the features at the same radius and 
different longitudes cannot be matched, implying that they are non-axisymmetric or time- 
variable or both. 

Horn & Cuzzi (1996) have measured the local power spectrum of the reflectivity of the 
B ring as a function of radius (from 92,000 km at its inner edge to 122,000 km at the Cassini 
division). In most regions the dominant wavelength varies between 100 and 200 km. In the 
outer 1500 km of the ring, and a few other isolated regions, there is substantial power at 
wavelengths as short as tens of km. There are also regions of 1000 km or so in which little 
or no irregular structure is visible. 

The particle-size distribution in Saturn’s rings is constrained by the Voyager occultation 
experiments (Zebker, Marouf, & Tyler 1985; Showalter & Nicholson 1990). These observa¬ 
tions suggest a broad distribution of particle sizes with an upper cutoff of a few meters, but 
cannot distinguish whether these particles are in differential rotation or locked by contact 
forces into a solid assembly. 


3. Theoretical models 


The irregular structure is surprising because viscous diffusion is expected to smooth out 
such structure on timescales much less than the age of the rings. The characteristic time 
required to smooth out structure on a radial scale Ar is 


(Ar)^ 4 

U = - -^ = 3x10* y 


Ar 


100 km 


100 cm^ sec * 


V 


( 1 ) 


where v is the kinematic viscosity, which can be estimated either from the rate of damping 
of density waves, or from kinetic theory together with estimates of the particle size and 
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optical depth. The very small value of ti, implies that all small-scale radial structure should 
be erased from the B ring, unless it is actively maintained by some mechanism. 

One possible explanation is that the irregular structure arises from gravitational shep¬ 
herding by small moonlets or large particles. However, objects large enough to influence the 
ring over scales > 100 km should also clear a gap around themselves, and such gaps are not 
seen (Esposito, O’Callaghan, & West 1983). 

One attractive hypothesis is that the irregular structure arises from local axisymmetric 
instabilities in the rings. There are two natural candidates for this instability: 


• Viscous instability. In most models for axisymmetric ring dynamics, the kinematic 
viscosity z/ is assumed to be a function only of the surface density jj,. The viscous 
instability then arises if the angular-momentum flux in the rings is a decreasing function 
of surface density, that is, if 


Ki = 


dn 


> 0 ; 


( 2 ) 


the growth rate of the instability is 7 = Sk^Ki, where k is the radial wavenumber. 
This is a secular instability, since the growth rate is proportional to the strength of the 
dissipative forces (i.e., 7 oc z/). If condition (2) is satisfied, ring material preferentially 
migrates from regions of low surface density to regions of high surface density, thereby 
adding mass to the high-density regions, and depleting the low-density regions even 
further. 


The viscous instability was hrst discussed in the context of accretion disks by Lightman 
& Eardley (1974), and was invoked to explain the irregular structure in Saturn’s rings 
by Lin & Bodenheimer (1981), Lukkari (1981), and Ward (1981). These authors were 
motivated by the observation that the viscosity in dilute rings appears to satisfy the 
instability condition (2). However, the B ring is unlikely to be dilute, unless the 
elasticity of its constituent particles is nnrealistically high. Both analytic kinetic theory 
and iV-body simulations of the viscosity in a dense ring composed of inelastic hard 
spheres show that the instability condition ( 2 ) is not satisfied by a wide margin (Araki 
& Tremaine 1986; Wisdom & Tremaine 1988); therefore, the viscous instability is 
nnlikely to operate in Saturn’s B ring. 


• Viscous overstability. Density waves can propagate in planetary rings, and dozens 
of density wavetrains have been identihed in Saturn’s A ring. The collective effect 
that dominates density-wave propagation is the self-gravity of the ring. Self-gravity 
is important, even though the ring is mnch less massive than the planet, because the 
thickness of the ring is also much less than its radius. A more precise statement is that 
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Toomre’s Q parameter is of order unity in the rings. Viscosity can lead to either decay 
or growth of density waves; the latter case, known as viscous overstability, occurs if 

K 2 = + |z/ - e > 0; (3) 

here u and ^ are the kinematic shear and bulk viscosity, and both are assumed to 
depend only on the surface density /i. The growth rate is 7 = Unfortunately, 

the assumptions that the stress tensor can be represented by a shear and bulk viscosity, 
and that these depend only on the surface density, are harder to justify than in the case 
of the viscous instability, since an overstable ring oscillates on a timescale comparable 
to the interparticle collision time. 

The existence of this instability was hrst pointed out by Kato (1978) for accretion 
disks (where the important collective effect is pressure rather than self-gravity), and 
by Goldreich & Tremaine (1978) for planetary rings. Borderies, Goldreich & Tremaine 
(1985) describe a crude kinetic theory for a ring composed of inelastic, closely packed, 
hard spheres, and show that the viscous instability is present in this approximation. 
Schmidt & Tscharnuter (1995) discuss in detail the local linear stability of a differen¬ 
tially rotating, two-dimensional, isothermal, fluid disk, in which the viscosity depends 
only on surface density, and the effects of self-gravity are included; they derive a cu¬ 
bic dispersion relation that contains the viscous instability and viscous overstability 
as special cases. Schmidt & Tscharnuter (1999) describe numerical calculations of 
the nonlinear evolution of these disks, which show that the viscous overstability can 
lead to a rich radial structure with surface-density contrast of order unity. Numer¬ 
ical experiments on self-gravitating collections of ring particles appear to exhibit a 
viscous overstability in some cases (Mosqneira 1996; Daisaka, Tanaka, & Ida 2001; 
Salo, Schmidt & Spahn 2001). However, the most unstable wavelengths in the linear 
calculations, and the wavelengths with signihcant power in the nonlinear calculations, 
are at most a few hundred meters—almost three orders of magnitude smaller than the 
dominant wavelengths in the irregular structure in the B ring (Horn & Guzzi 1996). 
Schmidt & Tscharnnter (1999) argue that nonlinear wave-wave interactions transfer 
the flnctnation power to much larger wavelengths, bnt so far this snggestion remains 
nntested. 

A quite different proposal is that some of the irregular structure arises from unstable 
ballistic transport (Lissauer 1984; Durisen et al. 1989, 1992; Durisen 1995). Impacts by in¬ 
terplanetary particles erode the ring particles and redistribute their mass to adjacent annnli. 
This process is unstable because high-density regions of the ring tend to absorb more of the 
ejecta than neighboring regions. However, the growth rate of the instability is negligible for 
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optical depths > 1.5. Thus, ballistic transport is probably unable to explain the irregular 
structure in the bulk of the B ring. 


4. Shear stress in dense planetary rings 

In this section we investigate the dynamics, stresses, and stability of solid and liquid 
ring phases. 


4.1. Equations of motion 


We examine a ring orbiting a point mass M. Test particles on circular orbits travel 
at the Keplerian angular speed ffii-(r) = (GM/r^Y^'^. We employ a rotating Cartesian 
coordinate system in which points radially outward, e^ points in the direction of rotation, 
and the origin, at radius R, rotates around the central mass at angular speed 12 = Qk{R)- 
We work at distances from the origin that are small compared to the orbital radius R (Hill’s 
approximation). We ignore the vertical structure of the disk, treating it as a razor-thin sheet 
with surface density /i(x, f). We restrict ourselves to axisymmetric disturbances, so that 
d/dy = 0. The Euler and continuity equations then read 


du du 
dt ^dx 


dv dv 
dt ^dx 


djj, d , . 


312^0; -|- 2Qv H--—, 

/i dx 


—212m -|- 


i dE^y 
fi dx ’ 


0 , 


( 4 ) 


where v(x, 2) = u{x,t)ex + v{x,t)ey is the velocity in the rotating frame, and = / (Jikdz 
is the vertically integrated stress tensor (dimensions of force per unit length). We define the 


tangential shear to be 


dv 

dx 


( 5 ) 


In a disk with zero stress gradients, a solution of the equations of motion is 

u{x,t) = 0, v{x,t) = vk{x) = —|12a;, (6) 

corresponding to circular Keplerian orbits and constant shear s = sk = ~|12- In a solid 
disk with zero shear, a solution of the equations of motion is 

dT.^x 


u{x,t) = 0, v{x,t) = Vs = constant, Ti^y = constant. 


dx 


y{2nvs + 3n^x). (7) 
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If we assume that the surface density is constant, that the zero-shear region extends from 
Xi to X 2 , and that the tensile stress vanishes at the edges of the solid region (i.e. = 0 at 

Xi,X 2 ) then the last of these equations can be integrated to yield 

- xi){x2-x), Vs =+ X 2 ). ( 8 ) 

Thus, any solid annulus is subject to a tensile stress, the maximum of which occurs at its 
midline and is equal to 

Sxx,max = (9) 

where Ax = X 2 — xi is the width of the annulus. 

If we linearize the first two of equations (4) with respect to a state of uniform shear 
{uq = 0, dvo/dx = s), and neglect perturbations in the stress tensor, then small disturbances 
are stable if and only if s > —2VL. This is the well-known Rayleigh criterion for the stability 
of Couette flow (Chandrasekhar 1961). 


4.2. Tensile stress 


Equation (9) for the maximum height-integrated tensile stress (force per unit length) in 
a solid annulus can be converted to an equation for the ordinary tensile stress (force per unit 
area) by approximating the vertical structure of the ring as that of a homogeneous slab with 
thickness h and density p = jj,/h (note that p is less than the density of the ring particles, 
by the hlling factor): 


a. 


-'tctc.max 


xx.mdix 


h 


= |pr2^(Ar)^ = 4x10^ dyn cm 


-2 


P 


0.3 g cm 


-3 


10 


10 


cm 


Ar 


100 km 


( 10 ) 

in which we have inserted parameters appropriate for Saturn. The maximum width of a solid 
annulus with a given tensile strength may be called its “tidal width”. If we identify the tidal 
width with half of the dominant wavelength of 100 km seen in the irregular structure (i.e. 
Ar = 50 km) then we require that a frozen assembly of ring particles have a tensile strength 
or yield stress > 1 x 10® dyn cm“^. 


The strength of ring-particle assemblies is, of course, very difficult to estimate. An 
upper limit is the yield stress of solid ice, which is Umax ~ 10^ dyn cm“^ at temperatures 
5-10 K below freezing (Schulson 1999), and probably higher at the much lower temperatures 
characteristic of Saturn’s rings. However, the actual strength of ice-particle assemblies can 
be many orders of magnitude lower: comets, which may be rubble piles of icy particles, 
have yield stresses Umax ^ 10® dyn cm“^, and Asphaug & Benz (1996) derive an even lower 
strength < 10^ dyn cm“^ for Comet Shoemaker-Levy 9. 



Experimental work on sticking of cm-sized ice particles is reviewed by Bridges et al. 
(1996). The experiments show no sticking unless the ice surfaces are coated with uncom¬ 
pacted frost; however, thin water-frost layers can lead to sticking forces up to 10^ dyn over 
contact areas ~ 0.01 cm^ (see also Supulver et al. 1997). On Saturn’s ring particles, frost is 
likely to be produced by dust-particle impacts, and removed by collisions with other particles; 
Hatzes et al. (1991) conclude that centimeter-sized or larger particles will be frost-coated, 
although this conclusion relies on a very uncertain estimate of the optical depth in grains. 
Whether or not frost layers are present on B-ring particles is therefore an open question; we 
shall assume that frost is present in order to have signihcant sticking forces. 

If the sticking force scales with the contact area, then the forces found by Bridges et al. 
(1996) would imply a tensile strength per unit area of 10® dyn cm“^. This is certainly an 
overestimate, since the contact area between ring particles is only a fraction of the total area. 
If there is a wide range of particle sizes, the contact points and the yield stress are likely to 
be dominated by small particles, which provide a kind of “cement” to hold the large particles 
together. We now make a crude estimate of the yield stress. For spherical particles of radius 
R, the contact area is A = 2TrRd where d is the thickness of the frost layer. We assume that 
the contact force is proportional to the contact area, that the force is ~ 10"^ dyn for a contact 
area A = 0.01 cm^, that the thickness of the frost layer is d ~ 10“^ cm, and that there are 
iV 1 cm^/Tri?^ contacts per cm^. Then the yield stress is Umax ~ 2 x 10^(1 cm/ R) dyn cm~^, 
comparable to the required yield stress if R ^ 0.2 cm. Note that most of the mass of the 
ring could be—and likely is—in much larger particles, while most of the sticking force comes 
from the small particles. 

Therefore, within the very large uncertainties, it is plausible that annular ring-particle 
assemblies as large as Ar ~ 50-100 km could have the required strength (umax ~ 10® dyn cm”' 
to survive tidal stresses. 

We now briefly discuss the orbital stability of solid annuli. Laplace and later Maxwell 
showed that a solid annulus orbiting a planet is unstable, with growth time 12/^2. How¬ 
ever, Fridman, Morozov, & Polyachenko (1984) point out that the Laplace-Maxwell analysis 
assumes that the annulus is completely rigid, in particular that the internal sound speed 
c ^ fir. In fact c fir — 20 km s”^ for any plausible material. In this case the nature of 
the instability changes dramatically (Fridman, Morozov, & Polyachenko 1984): the growth 
rate becomes much smaller, ^/3mc/r, where m is the azimuthal wavenumber; the instability 
is predominantly azimuthal; and the m = 1 instability in the rigid ring is replaced by an 
instability that is present for all m and grows faster as m increases. The physical basis for 
the instability is simple: if a mass element moves closer to the element in front of it, it 
is pushed backwards by elastic forces. Thus it loses angular momentum, its angular speed 
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increases, and it moves even closer to the element in front. 

One limitation of the analysis of Fridman, Morozov, & Polyachenko (1984) is the neglect 
of the rigidity of the ring: they assume that the modulus of compression is non-zero but the 
modulus of rigidity /i is zero. Rigidity suppresses the instability at short wavelengths, for 
m > merit = where p is the density in the ring. 

Thus there are several possibilities that may enable solid ring-particle assemblies to be 
stable. Possibly the rigidity is large enough to stabilize all azimuthal wavenumbers; possibly 
the solid regions are arcs, rather than annuli, spanning angles < 27r/mcrit; possibly there 
is a slowly growing instability that disrupts the solid annuli, so that they are constantly 
dissolving and re-forming; or perhaps other effects such as pressure from adjacent parts of 
the ring can stabilize a solid annulus. We cannot say which of these possibilities is more 
likely, but it would be premature to dismiss the possibility of solid annuli in the ring on the 
basis of stability arguments. 


4.3. Shear stress 


The simplest model for the shear stress is = Es, where E = f pdz, the vertically 
integrated viscosity, is assumed to be a function only of the surface density p. Many authors 
have investigated the stability and evolution of models of this kind (Lightman & Eardley 
1974; Goldreich & Tremaine 1978; Kato 1978; Lin & Bodenheimer 1981; Schmidt & Tschar- 
nuter 1995, 1999). However, granular systems generally behave as non-Newtonian fluids, in 
which the viscosity is a function of the shear, so that the shear stress llxy is likely to be a 
nonlinear function of the shear. 


For example, a simple kinetic theory for dense planetary rings is described by Borderies, 
Goldreich & Tremaine (1985). In their model, the shear stress and radial stress are 


^xy QiE 




sgn (s) 


V — 

^XX 


-q2F 




( 11 ) 


here qi and q 2 are dimensionless constants of order unity, and p is the density of the ring 
particles or the mean density inside the ring (the theory is not accurate enough to distinguish 
these two densities, as it assumes in effect that the ring is incompressible). The dimensionless 
factor 




P 


( 12 ) 


0.3 g cm“3/ Vl0^°cm/ 

is the enhancement in the vertical gravitational held due to the self-gravity of the ring. The 
negative radial stress Tixx represents the hydrostatic pressure required to maintain the ring 
thickness. 
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More generally, the shear stress is likely to be a complex function of the shear, the surface 
density, and perhaps also their histories. Analytic models for the constitutive relations of a 
dense ring of inelastically colliding particles are still in a primitive state, particularly when 
cohesive forces are present. On the other hand A^-body simulations of ring dynamics, which 
will be an indispensable guide to the correct analytic theory, have so far been restricted to 
the Keplerian shear rate and do not include the possibility of cohesion. 

Given our present ignorance, the only practical approach is to parametrize a fairly 
general set of constitutive relations, and investigate the behavior of the ring as a function of 
these parameters. We shall therefore assume that the shear stress and the radial stress are 
arbitrary nonlinear functions of the shear s and the surface density p. We can deduce some 
plausible general properties of these stresses. The shear stress s) should vanish if the 

shear vanishes (if there is no shear then there is no stress). We also expect that will 
become very large as the shear s approaches —2Q, since we have seen in §4.1 that circular 
orbits are unstable when s < —211. It is also natural to assume that has the same sign 
as the shear s; however, we do not assume that T,xy is monotonic in s, and in subsequent 
sections will focus on instabilities that can arise when Ti^y is a decreasing function of s in 
some interval. The radial stress s) will be negative if the ring behaves as a fluid that 

exerts a pressure force, but may become positive in regions of negligible shear when the ring 
material freezes. 


4.4. Linear stability 


As described above, we shall assume that the stresses are functions of the local shear 
and surface density, 

^xx f (h") "S) ) ^xy fl'(h') '®)- (13) 

We then linearize the equations of motion (4) around a state in which the surface density 
is uniform, the unperturbed motion is Keplerian (eq. 6), and the stresses are constant. We 
assume that the perturbations vary as exp{ikx + ■yt). We hnd 


yui — 2Vtvi 

yvi + \VLui 
7/ii + ikfioUi 


ik 

ho 

ik 


(//.hi + ikfsVi) , 
(gy^i + ikgsVi) , 


ho 

0 , 


( 14 ) 



-11 - 


where fs = {Of /ds)o, etc. The resulting dispersion relation is^ 

+ k^{fs9fM - fi^9s) - 2VLk'^ii9^ = 0. (15) 

We can recover the viscous instability and overstability of §3 by neglecting the radial 
stress (/ = 0) and assuming that the ring is a Newtonian fluid so that g = nsvi^jj). Equation 
(15) then simplihes to 

^3 + ^ ^ = 0 . (16) 

dfi 

If the viscosity is small, this equation has two possible solutions: either 7 = —3k‘^d{g,h')/dg,+ 
0 (z/^), corresponding to the viscous instability (eq. 2 ), or 7 = iifl + k‘^[^i^{diy/dfi) + u], 
corresponding to the viscous overstability (eq. 3, with the bulk viscosity ^ so that 

the radial stress vanishes). 

Thus the viscous instability and overstability arise in the limit where the stresses f,9^ 
0 at hxed wavenumber k. In this paper, we will focus instead on the short-wavelength limit, 
in which k ^ 00 at hxed values of /, g. In this case the roots of the dispersion relation are 
given by 

J ^-fsi + 0(k«), T^ = p(/„-/.^) +0(*:«). (17) 

\ 9sJ 

The hrst of these is the largest in absolute value, and leads to a rapidly growing instability 
if and only if gs < 0 . This instability does not arise in a Newtonian disk, in which g = Es 
so that gs = E > t). 

Assuming that the stresses / and g are comparable in magnitude, and that /* f/kl 
and ///i, with similar relations for g, the approximations that lead to the formula (17) 
are valid so long as k'^g > If we generalize equation (11) to 

g{s,f^) = E^w{s), (18) 

where w'{s) is of order unity, then this condition becomes A < 27rF^/^/i, where h = p/p 
is the effective thickness, and we have replaced the wavenumber k by For typical 

values of E (eq. 12), the condition becomes A < 20h. We expect that our two-dimensional 
approximation remains valid so long as A > h, so there is roughly one decade in wavenumber 
in which the instability should be present whenever w'{sk) < 0 . 


^This analysis neglects the effects of the ring self-gravity; these can easily be added to (15) and have no 
important effect on our conclusions. 
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We have established that a strong short-wavelength instability is present in rings in 
which dg/ds < 0, that is, in which the shear stress becomes larger in absolute value when 
the shear rate becomes smaller in absolute value. We now analyze the nonlinear behavior of 
this instability in a ring model that is simplihed even further. 


5. A toy model 


In this section we investigate a simple fluid system that illustrates an instability similar 
to the one we derived in the preceding section. Consider a homogeneous two-dimensional 
system of viscous incompressible fluid with surface density /i. The system is inhnite in the 
y-direction, and all of its physical properties are assumed to be independent of y. The 
system is in shear flow in the y-direction, described by the velocity field v = (0, v{x, t)). The 
system has periodic shearing boundary conditions in the x-direction; that is, v{x + Ax, t) = 
v{x, t)+SKAx, where Ax is the period in x, sk is the mean shear (i.e. the shear averaged over 
length scales much larger than Ax). The only non-trivial component of the Euler equation 
reads 

dv 1 da^y 
dt y dx ' 

where a is the stress tensor. 



We let s{x,t) = dv/dx denote the shear. We shall assume that the stress is a nonlinear 
function of the shear, a^y = g{s). Differentiating equation (19) with respect to x, we obtain 


ds 

m 


1 

fx dx‘^ 


gi.s). 


The periodic boundary conditions require that 


( 20 ) 


r*xH-Aa: 


s{x, t)dx = skAx for all x. 


( 21 ) 


Initially we assume that the system is in a state of uniform shear, s{x,t) = sk- A 
linearized analysis of equation (20) reveals that small disturbances of the form exp( 7 t-|-i/i;a;) 
satisfy the dispersion relation 

^ k^g'jsK) 

fi 

exactly the same as the largest root in equation (17). Thus we have reproduced the instability 
derived in §4.4 in a simpler system, in which we can more easily investigate the nonlinear 
evolution of the instability. 
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The final state of an unstable system with g'{sK) < 0 is straightforward to describe. 
First we define 


G(5) ^ / g{s')ds\ Fit) ^ 


r*ai+Aai 


G(s(a;, t))dx. 


Using equation (20) and the boundary conditions, we have 


(23) 


dt 


f-x+Ax o 1 

g{s)—dx = - 
j. at /i 


rx+Ax q2 

9 {s)^gis) dx 




2 

dx < 0. 


(24) 


Since F decreases with time, the system must approach a steady state as t > cx). The 
steady-state solutions of equation (20) are those in which gis) is a linear function of x] the 
periodic boundary conditions require that the linear term vanishes, so g{s) = constant in the 
steady state. However, the only solution of this form that satisfies the mean shear constraint 
(21) is s = sk, which is unstable. The resolution of this apparent contradiction is that the 
final state of the unstable system is piecewise constant in the shear s, alternating between 
states s_ and that satisfy 5 '(s_) = gis+) = gisx)- The fractions f± of the spatial interval 
Ax that are occupied by each of these two phases are constrained by the relations 


/- + /+ = 1, f-s- + f+s+ = Sk, (25) 

where the second equality follows from equation (21). However, the equations provide no 
information on the spatial distribution of the two phases in the final state; for example, the 
characteristic domain size is unknown. 

The evolution from a uniform but unstable initial state into two distinct phases is known 
as spinodal decomposition in studies of alloys and binary fluids; in this case the shear s is the 
concentration of one of the two components, g{s) is the chemical potential, the functional 
F is the free energy of the system, and equation (21) corresponds to conservation of mass 
(Bray 1994). This conservation law places our model in the class of phase transitions with 
a conserved order parameter. 

The function gis) must satisfy several constraints in any physically plausible system; 
(i) gis) is an odd function of s (since the directions +y and —y are equivalent); (ii) 5 ^( 0 ) = 0 
(since there is no stress if there is no shear); (iii) gis) has the same sign as s (since the 
viscosity is positive). A functional form that is general enough for our purposes is 

5f(s) = s(as^-|-&|s|-|-c); (26) 

condition (iii) then requires that a > 0, c > 0, and b > — (4ac)^/^. 
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With this form for g{s), there is an unstable region {g'{s) < 0) if and only if 6 < 
— (3ac)^/^. When this condition is satished, all values of |s| between > 0 and Sb > 0 are 
unstable, where 

b (P — 3ac)^/^ 

Sa,b 7, ^ 7, 

3a 3a 

The interval Sa < |s| < s;, in which a uniform shear is unstable is called the spinodal interval 
(Figure 1). 

There is also interesting behavior outside the spinodal interval. Suppose that the system 
is in uniform shear Si outside the spinodal interval, Si > Sb- Suppose that there is another 
shear state outside the spinodal interval, S 2 < Sa, such that g{si) = g{s 2 ). Then if a small 
element 6x changes its shear state from si to S 2 , and the remaining fluid increases its shear 
rate from Si to Si + (5si so as to satisfy the mean shear constraint (21), the resulting change 
in free energy is 

SF = (5a;[G(s2) - G(si) - G"(si)(s2 - si)], (28) 

which is negative if 

G'K) = sK) > , (29) 

-52 ~ -Si 

Shear states in which this inequality is satished—the interval Sp < s < Sa and s;, < s < Sg in 
Figure 1—are metastable, because the mixture of two distinct phases has a lower free energy. 

Equation (20) is ill-posed, since the growth rate of the instability (22) becomes extremely 
large for short-wavelength disturbances. It is convenient, and physically plausible, to mitigate 
this violent instability by modifying equation (20) to 

ds 1 (9^ d'^s 

dt ^ ^dx'^’’ 

where e is a small positive constant related to the thickness of the interface. The additional 
term can be thought of as penalizing the growth of solutions with large gradients. 

We now investigate the properties of solutions of equation (30). 

It is simple to see that the mean shear is conserved, 

x-\-/lx 

s{x,t)dx = 0, (31) 

which implies that if the mean shear condition (21) is satished initially then it is automatically 
satished for all time. 
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Fig. 1.— (left panel) A possible form for the shear stress g{s) as a function of shear s, in the 
model of §5. Uniform shear is unstable in the interval Sa < s < Sb (eq. 22). (right panel) The 
function G{s) dehned in equation (23), shown for the same parameters and on the same scale 
as the left panel. States in the interval Sp < s < Sa and Sb < s < Sg are metastable. The 
values Sp and Sq are determined by solving the equality corresponding to the inequality (29), 
along with the condition g{sp) = g{sq). The solution of these two simultaneous equations is 
equivalent to hnding the straight line tangent to G(<s) at two points. 
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thus only wavenumbers that satisfy are unstable; the most unstable 

wavenumber is given by 


Let us dehne 


f*X-\-Ax 


F(t) = 


G{s{x,t)) ^ (ds 

h 2 ^ 


/i 


dx 


dx, 


(33) 


where G(s) is defined in equation (23). It is straightforward to show for periodic boundary 
conditions that 


dF 

dt 


(^x-\-Ax 


fi dx 


d^s 

dx^ 


dx < 0. 


(34) 


Since F decreases with time, the system must approach a steady state a.s t ^ oo. 


The steady-state solutions of equation (30) are those in which g{s)/g, — e{d‘^s/dx‘^) is a 
linear function of x; the periodic boundary conditions require that the linear term vanishes, 
so the steady state solutions satisfy 


dx'^ 


1 , , 

—g{s)+p = 0, 
e/i 


(35) 


where p is a constant. This equation is equivalent to motion in the potential —G(s)/(e/i) +ps, 
so there is an “energy” integral 


1 

2 



— G(s) +ps = q; 
e/i 


(36) 


the properties of the steady-state solutions can thus be determined by examining the contours 
of the left-hand side of (36) in the phase plane with coordinates {s,ds/dx). 


We now introduce dimensionless variables. 






1/2 


epj 


X, 


r = 

e/i^ 


y = 


a\i/2 


- 1 

C/ 


(37) 


recall that a and c are positive. Equation (30), with the stress-shear relation (26), becomes 

= + B (38) 

The boundary conditions are periodic, with period A,^ = (c/e/i)^/^Aa;. Equation (38) is the 
Cahn-Hilliard equation, which has been widely used to model pattern formation in phase 
transitions. 


Consider small disturbances to uniform shear flow with shear sk- The dispersion relation 
for disturbances of the form exp( 7 r -|- in^) is 


7 = -K^(l + 2B|j/ji-| + 3j//) - ft', 


( 39 ) 



(40) 


where yx = {cl/cY^'^sk- There is instability if 7 > 0 for some k, which occurs if 


B < 


1 + 3 |/^ 

2 ||/ic| 


If this instability is present, the stress as a function of shear g{s) is triple-valued: there are 
three roots < sk < to the equation g{s) = g^sx)- The root sx has g\sx) < 0 and 
hence is unstable, while g'{s±) > 0 so the roots are stable. 

Almost all of the steady-state solutions of equation (38) are bounded and periodic, and 
can be expressed analytically in terms of Jacobian elliptic functions Novick-Cohen & Segel 
(1984). However, these are not the final state of the system: it turns out that all of the 
periodic solutions are unstable (Carr Curtin & Slemrod 1984). The only stable, stationary 
solution is the “kink” solution. 


y{x) = ± ~ l)^'^^tanh 


(lr>2 i\l/2 X 


(41) 


where the two ± signs are independent and we require that B < —3/v^- Of course, this 
solution does not satisfy our periodic boundary conditions. 

Numerical integration of the partial differential equation (38) for a unstable initial state 
shows that the system evolves to a state in which the shear is almost always nearly equal 
to either s_ or s+, just like the solutions of the simpler equation (20). Moreover, the 
interfaces between high- and low-shear domains gradually drift, so that high-shear and low- 
shear domains eventually coalesce. Once the distance between interfaces is large compared 
to unity, the shape of the shear curve in the transition region is given approximately by 
the kink solution (41). Because of the conservation law (31), the interfaces cannot move 
independently. Their interactions can be modeled by treating each interface as a particle 
that exerts forces on other particles (Kawasaki & Ohta 1982; Majumdar & Huse 1995). This 
process is known as coarsening or ripening. In an inhnite system, the coarsening process 
continues indehnitely, although at a slower and slower rate as the domain sizes grow. Thus 
complete thermodynamic equilibrium is never achieved. 

The coarsening process in the one-dimensional Cahn-Hilliard equation exhibits scaling 
behavior, that is, the characteristic distance between the interfaces grows indehnitely but 
the correlation function retains the same shape. However, in a planetary ring the coarsening 
will be halted when the distance between interfaces becomes comparable to the tidal width 
of §4.2. Thus our toy model leads naturally to the conclusion that the irregular structure 
in Saturn’s rings should have a characteristic size comparable to the tidal width, which we 
estimated could be as large as ~ 100 km. 
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Assuming that the characteristic width of the interface between solid and liquid phases 
is L, the parameter e is of order L‘^u where u is the kinematic viscosity. It is plausible that 
L is comparable to the radius of large ring particles. 

One interesting unresolved issue is whether our one-dimensional approximation is ade¬ 
quate. Coarsening in the Cahn-Hilliard equation is much faster for systems with more than 
one dimension, where the free energy associated with the interfaces or domain walls can be 
reduced by reducing their radius of curvature. Is the modest curvature of annuli in the B-ring 
sufficient to make this an important contributor to coarsening dynamics? A second unre¬ 
solved issue is the importance of noise. The Cahn-Hilliard equation is a zero-temperature 
model; adding thermal noise converts the equation to a hnite-temperature model. Is such a 
model more relevant for Saturn’s B ring, and if so, what sets the effective temperature (im¬ 
pacts of small bodies? gravitational wakes from large ring particles?)? Are the fluctuations 
sufficient to trigger phase separation in the metastable region of the Cahn-Hilliard equation? 


6. Discussion 

We have suggested that the irregular structure in Saturn’s B ring arises from the forma¬ 
tion of solid ring-particle assemblies, which are limited in size by the competition between 
tidal forces and the tensile strength of these assemblies. We have shown that if the shear 
stress is a decreasing function of the shear, dTi^y/ds < 0, then the ring is unstable at short 
wavelengths. Using a toy model of an incompressible, non-Newtonian fluid in shear flow, we 
have argued that this instability leads to the formation of domains in which the shear takes 
on one of two values; it is natural to identify one of these with the solid phase, and one with 
a liquid phase in which the (absolute value of the) shear rate is greater than Keplerian. The 
toy model is ill-posed at short wavelengths; if this is remedied by adding a term that can 
be associated with the energy of domain walls (the term proportional to e in eq. 30), then 
the domains gradually coarsen, leading to structure on scales much larger than the original 
unstable wavelengths. 

There are many shortcomings in this work. We have not derived a constitutive rela¬ 
tion for a ring-particle fluid with cohesion, and thus we have not shown that the instability 
condition dJ^^y/ds < 0 is satished in rings with cohesive forces (although the simple kinetic- 
theory model of Borderies, Goldreich & Tremaine (1985) is marginally unstable; see eq. 11). 
We have not derived the statistical properties of the equilibrium ring structure, in which 
the growth of larger domains due to coarsening is balanced by their destruction from tidal 
forces. We have not discussed why the rigidly rotating solid regions should have different 
phase functions or albedo from the regions in which the particles are in differential rotation. 
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We have not shown that solid annuli are dynamically stable. Finally, we have not estab¬ 
lished that ring particles stick to one another with the required cohesive strength; laboratory 
experiments suggest that a frost layer is needed but we do not know whether such layers are 
present in the B ring. 

The most powerful tool we have to investigate at least some of these issues is iV-body 
simulations of local ring dynamics. Such simulations have already provided considerable 
insight into the dynamics of rings composed of inelastically colliding particles (Wisdom & 
Tremaine 1988; Richardson 1994; Mosqueira 1996; Salo, Schmidt & Spahn 2001); what is 
required is to generalize them to include cohesion. 


I thank David Huse and Roman Rafikov for helpful discussions, and Ignacio Mosqueira 
for detailed and thoughtful comments. This research was supported in part by NASA grant 


NAG5-10456. 
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